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We focus on Hofava-Lifshitz (HL) theory of gravity, and, in particular, on the Kehagias 
and Sfetsos's solution that is the analog of Schwarzschild black hole of General Relativity. 
In the weak-field and slow-motion approximation we analytically work out the secular 
precession of the longitude of the pericentre of a test particle induced by this solution. 
Its analytical form is different from that of the general relativistic Einstein's pericentre 
precession. Then, we compare it to the latest determinations of the corrections Aro to 
the standard Newtonian/Einsteinian planetary perihelion precessions recently estimated 
by E.V. Pitjeva with the EPM2008 ephemerides. It turns out that the planets of the solar 
system, taken singularly one at a time, allow to put lower bounds on the adimensional 
HL parameter i/iq of the order of 10"^'^ (Mercury) — 10"^'* (Pluto). They are not able to 
account for the Pioneer anomalous acceleration for r > 20 AU. 
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1. Introduction 

Hof avgPZI Jias recently proposed a four-dimensional theory of gravity that can be 
thought of as a candidate for a UV completion of General Relativity (GR), since it 
is power-counting renormalizable. This theory admits the Lifshitz scale invariance: 
X — 6x, t b'^t, and, after this, it is referred to as Hofava-Lifshitz (HL) theory. 
It is g = 3 in the UV, and hence HL theory exhibits anisotropy between space and 
time, while relativistic scaling with g = 1 is recovered in the IR regime. Actually, in 
the original formulation of HL gravity there are problems in recovering GR at large 
distances, unless the condition of "detailed balance" (used by Hofava to restrict 
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the number of po ssible parameters) is relaxeci^El. Thus, phenomenologicahy viable 
modificationa2lSI of th e the ory have been considered, even though some problems 
still seem to be presenlPI^l 

Since its formulation, HL gravity has been deeply investigated; several of its 
applications (e.g. cosmology, study of exact solutions, lensing, black holes and ther- 
modynamics of black holes) have been analyzed (see Refs. [U [10] and references 
therein). In particular, a static, spherically symmetric solution has been founc3by 
Kchagias and Sfetsos (KS). Such a solution is the analog of Schwarzschild one of 
GR and, moreover, it asymptotically reproduces the usual behavior of Schwarzschild 
spacetime. The KS solution is obtained without requiring the projectability condi- 
tion, assumed in the original HL theory. Sph erically sy mmetric solutions with the 
projectability condition are, however, aval It is important to point out 

that, up to now, in HL gravity the gravitational field is purely geometrical: in other 
words, the way matter has to be embedded still needs to be studied. Furthermore, 
since the way matter couples to gravity is not clear, the fact that particles move 
along geodesic is not granted, and deviations from geodesic motion are expectecff^l 
That being said, it is possible to consider the KS solution as a toy model which 
can be studied to obtain a better theoretical insight into the fundamental aspects 
of HL gravity because of its relative simplicity. Actually, this has already been 
done, focusing on ir npli cations of this solution in various physical scenarios, ranging 
from accretion diskal^ to black holes thermodynamics (see Ref. [ini and references 
therein). 

Since the Newtonian and lowest order post-Newtonian limits of the KS solution 
coincides with those of GR, it is manifest that this solution of HL gravity is in 
agreement with the classical tests of GR, as it has recently been shown in Ref. 17 In 
this paper we focus on the first corrections to the GR behavior, in order to evaluate 
the impact of these corrections on solar system orbital motions. Furthermore, we use 
data from solar system observations to constrain the KS solution. In particular, we 
want to investigate the weak- field and slow-motion approximation of the KS solution 
by working out the induced pericenter precession of a test particle (Section [2]) and 
comparing it to recent observations-based determinations of the non-standard rates 
of the perihelia of some planets of the solar system (Section |3|). Section |4] is devoted 
to the discussion and the conclusions. 

2. Calculating the perihelion precession 

In order to study orbital motion in the gravitational field of the Sun, we start from 
a static and spherically symmetric metric 

with = d-d"^ + sin ddip^ , in Schwarzschild- like coordinates. In particular, for the 
asymptotically flat KS solution we havJ^l 



18, 2010 2:1 WSPC/INSTRUCTION FILE HLIJMPA 



Phenomenological Constraints on Hofava-Lifshitz Solution 3 

with M = ^,so that 

(3) 

Let us re-write e"^"^^ as 



e''^'') = 1 + il^r 



1_ l + — 



(4) 



Let us focus on a typical Sun-planet system so that r « 1 AU; it is clear that some 
limiting condition on ip must hold so that 

r 

as in the Schwarzschild case. After posing ip = it turns out that it must be 

Vo»4(^^^ «4x 10-24 (6) 
for M = Mq and r « 1 AU. Indeed, in this case, one has 

e^^^)«l-^ + ^f:^V, (7) 



2M 2M^^ 



ds^ = (I- — + -j-^] c^dt^ - — ^—^dr'^-r'^d^'^. (8) 



so that 



Even though here we are concerned on the KS solution of HL gravity, our approach 
apply as well to an arbitrary perturbation of the Schwarzschild solution in the form 

e''W-l-^ + ^, (9) 

where a is a suitable perturbation parameter: in fact, on setting t/jq — A^^/a, it is 
possible to get constraints on a starting from the ones obtained on -00, knowing the 
mass of the source M = c^M./G. 

By inspection of the metric ([8]), it is clear that the first ^/jq— type corrections to 
the GR behaviour for large distances are of order r~^\ We aim at investigating the 
effects of these corrections that depend on terms in the form and we neglect 
the higher order effects of Schwarzschild metric (i.e. the terms that are proportional 
to (M)\fc>2). 



To this end, we introduce the isotropic radial coordinate r 

1 



r^ril + (10) 
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and substitute in ([S]). After approximating, the metric becomes 

ds^ =1 \ c^dr - iH 7 (dr^ + r^d^^) (11) 

or, introducing Cartesian coordinates such that = + + z^, 

^2 A 2X 2Xn / 2X 1X*\.^2 J 2 J 2N 

ds - 1 1- + c^dr - 1 + ~ t:—^ [dx^ + dir + dz^) ; 12) 

\ r woi' J \ r 2 ipor* ) 

in the following we will re-name f as rQ 

From eq. (fT2l) and from the post-Newtonian equations of motion of a test 
particlePl 

a;* = --^c^hoo^i - -c^hikhoo^k + hoo,ki''i' + {htk,m - -jhkrn,i^ i'^i™, i = 1,2,3, 

(13) 

valid for time-independent space-time metrics without off-diagonal terms, it is pos- 
sible to obtain the following expression for the ■i/;o~dependent acceleration 

I Zl -U - 

„2 



''''' Voc^rS 
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(14) 



where f = ^ is the planet's unit vector in the radial direction, v is the planet's 
velocity, and Vr is its component in the radial direction, i.e. Vr = v ■ f. Since typical 
solar systems's planetary speeds are of the order oi v ~ na ~ 10'* m s~^, where a is 
the semi-major axis and n = GM /a^ is the Keplerian mean motion related to the 
orbital period by Pb == 27r/n, we can neglect the terms and, especially, (^) w, 
which is even smaller since Vr is proportional to the planet's eccentricity ejSI. Thus, 
A^ji reduces to the entirely radial term 

4G4m4 ^ 

whose magnitude, in the case M = Mq and r w 1 AU, is quite smaller than the 
Newtonian monopole for t/jq > 4 (^) ; as we will see later, such a condition is fully 
satisfied in the solar system. As a consequence, it is possible to apply the standard 
Gauss perturbative approaclPQl to derive the secular, i.e. averaged over one orbital 



period, precession of the longitude of the periheliorPS w = + w, where $7 is the 
longitude of the ascending node and w is the argument of the perihelion. The Gauss 
equation for the node's variation ii^ 

tl^ . A„(-)sin(c. + /), (16) 

nay 1 — & sm7 ^a/ 



^It useful to remember that the isotropic form of Schwarzschild metric is, up to 

,2 / 2M 2M^\ 2,2 3A42\ . 2 2 ,^2n 



, 2 
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where An is the component of the perturbing acceleration A, whatever its physical 
origin may be, normal to the orbital plane, / is the inclinatior0 to the ecliptic, and 
/ is the true anomaly reckoning the planet's position from the the periheliorP^. 
The Gauss equation for the variation of uj iJ^Ol 



nae 



-ArCosJ + Ar[l + -] sin / 



cosm, (17) 



where p = a(l — e^) is the semi-latus rectum, and A^ and At are the components 
of A along the radial and transverse directions, respectively, both in-plane. The 
un-perturbed Keplerian elli pse , on which the right-hand-sides of eq. (fTB]) and eq. 
p?)) have to be evaluated, 

(18) 



1 + e cos / ' 

while the average over one orbital period has to be performed by means oil^ 

dt = ^ ^ (19) 
n(l -I- ecos jY 

Since, in this case, there is no normal component of A^^, the node is left unaffected. 
The calculations with eq. pT]) are made simpler by the absence of the transverse 
component of A^^ as well. As a consequence, we obtain that the advance of the 
perihelion per orbit is 

= -^on2a6(l-e2)3' ^^0) 
thus, the secular precession of tu, reintroducing physical units, is 

i^) -^A::^r"U ^^GM)\A + ^) ^ 3(GM)V2(4 -f e^) 

^ ''1'° \2ti) 2V'oc6na6(l-e2)3 2il;QC^a^/^{l - e'^f ^ ' 

Note that eq. (|21l) is analytically different from the gravitoelectric 0{c'^) Einstein 
precession which is proportional to {GMY/'^a~^/'^{l — e^)"^. 



2.1. Alternative derivation of the perihelion precession 

Here we will derive the perihelion precession due to HL gravity with an alternative 
approach based on the Lagrange planetary equationi^Ql valid when the perturbation 
considered can be derived from a potential function. Firstly, it is required to compute 
the average over one orbital revolution of the correction A[/ to the Newtonian 
potential U = —GM/r responsible of the perturbation. 



^Q,,ui,I can be thought as the three Euler angles determining the spatial orientation of a rigid 
body, i.e. the un-perturbed Keplerian ellipse which changes neither its size, fixed by a, nor its 
shape, fixed by e, with respect to a locally quasi-inertial frame. 



18, 2010 2:1 WSPC/INSTRUCTION FILE HLIJMPA 



= (22) 



6 L. lorio & M. L. Ruggiero 

In our case, from the expression of /iqo in eq. (jl2p it follows 

It is straightforward to obtain 

^ _ (GM)4(2 + e^) 

The Lagrange equations of the variation of uj and Q due to a generic perturbation 

are 

'1 -e2\ dn cosldn] 



1 



na 



de sin/ 9/ 



(24) 



= / . gg. (25) 

na^V 1 — sin/ «/ 



Since TZ^po does not contain the inclination /, we simply have 



which yields just eq. ((2T|) . 



3. Confrontation with the observations 



The expression of eq. (pT|) can be compared to the latest determinations of the cor- 
rections Aw to the standard Newtonian/ Einsteinian periehlion precessions of some 
planets of the solar system recently estimated by contrasting large planetary data 
sets of various types with accurate dyna mical forc e models of (almost) all known 
Newtonian and general relativistic In view of the relatively short 

orbital periods of some of the planets considered with respect to the time interval 
spanned by the observational data set used, i.e. about one century, the estimated 
Ari7 can be meaningfully compared to the predicted secular, i.e. averaged over one 
orbital revolution, precessions. Latest results obtained by PitjevgP^I a,re shown in 
Table [Hand Table [21 they are based on the EPM2008 ephemerides constructed from 
a data set of 550000 observations ranging from 1913 to 2008. It is interesting to note 
that she estimated corrections to the perihelion precessions of Neptune and Pluto 
as well, although the available observations do not yet cover full orbital revolutions 
for them. 

Since the corrections to the standard perihelion precessions of Table [T]- Table 
[2| are statistically compatible with zero, attributing them to the action of the HL 
acceleration of eq. ()14p implies that the maximum value of tJjq is infinity. Lower 
limits to it can be obtained from 

3(GM)7/2(4 + e2) 



(27) 



2[Aw ± (5(ATi7)],naxC6a9/2(l - e^)^ 

They are shown in Table [3] for the inner planets and Table [4| for the outer planets 
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The larger value is 8 x 10~^^ and comes from Mercury. As an a-posteriori check of 
the consistency of our perturbative results with the assumption of eq. ©, we note 
that our lower bounds on V'o are always larger than A{M/rY for all the planets 
considered. 

Let us note that -00 = 8 x 10~^^ (see Ref. [I7| yields for Mercury a precession of 
—0.07007 arcsec cty~^, which is quite different from the standard Einstein precession 
of +42.98 arcsec cty~^. 

Finally, we mention that since A^^ is a radial acceleration, although directed 
outward the Sun and not spatially uniform, it may be interesting to see if our 
limiting values for ipQ may, at least, explain the magnitude of the anomalous Pioneer 
acceleration (Apjo = 8.74 x 10"^" m s^"'^) experienced by the two twin probes after 
20 M^. The answer is negative. Indeed, the upper bounds on corresponding 
to the lower bounds on V'o obtained from Uranus, Neptune and Pluto orbiting in 
the spatial regions in which the Pioneer anomaly manifested itself in its presently 
known form, are smaller by about one order of magnitude. They are shown in Table 

El 

4. Discussion and conclusions 

In this paper we considered the spherically symmetric KS solution, which is a solu- 
tion of HL gravity, obtained relaxing both the detailed balance and projectability 
conditions. The Newtonian and lowest order post-Newtonian limits of KS solution 
coincide with those of GR, and we focused on the first corrections to the GR be- 
haviour, in order to evaluate the impact of these corrections on the solar system 
orbital motions and use the available data to constrain this particular solution of 
HL gravity. Actually, our results apply as well to arbitrary perturbations of the 
Schwarzschild proportional to l/r*. 

After having analytically calculated the pericentre precession of a test particle 
averaged over one orbital revolution with both the Gauss and Lagrange perturbative 
schemes, we compared our theoretical prediction to the corrections to the standard 
Newtonian/Einsteinian perihelion precessions of all the planets of the solar system, 
recently estimated by E.V. Pitjeva by fitting the dynamical models of the EPM2008 
ephemerides to 550000 observations of several types spanning the time interval 1913- 
2008. The upper bounds on the estimated corrections to the perihelion precessions of 
all the planets, taken singularly one at a time, yield lower bounds on -00 ■ The largest 
one occurs for Mercury, being of the order of 10""'^^, while Pluto yields the smallest 
one which is about 10"^**. Concerning Mercury, the corresponding lower bound for 
■00 yields a retrograde perihelion precession of at most —0.07007 arcsec cty~^, quite 
different from the standard general relativistic one. The lower bounds on ip^ from 
Uranus, Neptune and Pluto, moving where the Pioneer anomaly manifested itself in 
its presently known form, yield upper bounds on the corresponding HL acceleration 
one order of magnitude smaller than the anomalous Pioneer acceleration of 8.74 x 
10"i° m s-2. 
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If and when other teams of astronomers wiU estimate their own corrections to the 
standard rates of the planetary of periheha, it will be possible to repeat the present 
analysis. A further, complementary approach that could be followed, although very 
time-consuming, consists in re-processing all the planetary data sets with the fully 
modeled KS metric and looking at the values of the estimated parameters, including 
also tpo itself. 
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Table 1. Estimated corrections Aro, in arcsec cty~^ , to the standard 
Newtonian/Einsteinian secular, i.e. averaged over one orbital revolu- 
tion, precessions of the longitude of the p erih elia of the inner planets 
obtained with the EPM2008 ephemerideJ^SI, The quoted errors are 
not the mere formal, statistical ones. The corrections Aro account for 
any un-modelled dynamical effects, classical (i.e. Newtonian/general 
relativistic) or not. 



Mercury 


Venus 


Earth 


Mars 


-0.0040 ± 0.0050 


0.024 ± 0.033 


0.006 ± 0.007 


-0.007 ±0.007 



Table 2. Estimated corrections Aro, in arcsec cty"-*^, to the standard New- 
tonian/Einsteinian secular precessions of the longitude of the perihelia of the 
outer planets obtained with the EPM2008 ephemerideJSSl The quoted errors 
are not the mere formal, statistical ones. The corrections Act account for any 
un-modelled dynamical effects, classical (i.e. Newtonian/general relativistic) 
or not. 



Jupiter 


Saturn 


Uranus 


Neptune 


Pluto 


0.067 ±0.093 


-0.010 ±0.015 


-3.89 ±3.90 


-4.44 ±5.40 


2.84 ±4.51 



Table 3. Lower bounds on \tpo \ obtained from the 
corrections A-cu to the standard perihelion preces- 
sions of the inner planets of Table [T] by using the 
predicted rate of eq. II21I I. 



Mercury 


Venus 


Earth 


Mars 


8 X 10-^2 


6 X 10"" 


6 X 10"" 


8 X 10-15 



Table 4. Lower bounds on \ipo\ obtained from the corrections 
Aro to the standard perihelion precessions of the outer planets 
of Table [2] by using the predicted rate of eq. I I21II . 



Jupiter 


Saturn 


Uranus 


Neptune 


Pluto 


3 X 10-18 


1 X 10-1** 


2 X 10-22 


2 X 10-23 


8 X 10-24 



Table 5. Upper bounds of the HL acceleration of eq. II14I I 
at heliocentric distances of Uranus, Neptune and Pluto com- 
puted with the lower bounds on i/)o according to Table |4] 



20 AU 


30 AU 


40 AU 


5.4 X 10-11 m s-2 


5.5 X 10-11 m s-2 


3 X 10-11 m s-2 



